program chem_freestream_test
   use mod_chem
   use iso_fortran_env, only: dp => real64
   implicit none

   type(type_chem_phases) :: phases
   integer, parameter :: n = 1981
   real(dp), allocatable :: T(:), X(:, :), rho_compared(:), H, U
   real(dp) :: p = 101325.0_dp, rho

   integer :: l, i, ll

   integer :: file_io(100)

   phases = chem_file_load("air_5species_NASA9.yaml", "air")

   print*, eos%gas%get_mass_fractions()


end